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ABSTRACT 

The thermal stability of a weakly magnetized, rotating, stratified, optically thin plasma 
is studied by means of linear-perturbation analysis. We derive dispersion relations 
and criteria for stability against axisymmetric perturbations that generalize previous 
results on either non-rotating or unmagnetized fluids. The implications for the hot 
atmospheres of galaxies and galaxy clusters are discussed. 
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C3 1 1 INTRODUCTION 

Galaxies and galaxy clusters are embedded in hot atmospheres of virial-temperature gas. The evolution of these systems 

■ depends crucially on whether these atmospheres are subject to thermal instability (Parker 1953; Field 1965; Defouw 1970), 
^ , so that condensations of cold gas can form and grow as a consequence of radiative cooling. Motivated by this astrophysical 

■ question, much effort has been devoted to studying the thermal stability of a stratified, optically thin plasma by means of both 
00 , linear-perturbation analysis (Malagoli, Rosner, & Bodo 1987; Loewenstein 1990; Balbus 1991; Binney, Nipoti & Fraternali 
C*") ■ 2009; Balbus & Reynolds 2008, 2010; Nipoti 2010, hereafter N10) and numerical hydrodynamical simulations (Kaufmann et 
\Q \ al. 2009; Joung, Bryan, & Putman 2012; McCourt et al. 2012; Sharma et al. 2012). The thermal-stability properties of the 
C^D ■ fluid are influenced by different physical processes. For instance, unmagnetized media behave differently from even weakly 

magnetized media. Similarly, the stability properties depend significantly on whether the fluid rotates and on the specific 
rotation law. While the effect of magnetic fields is accounted for in several of the aforementioned works, less attention has 
been paid to the thermal stability of a rotating stratified plasma. The effect of rotation on the thermal stability, which was 
discussed in simple cases by Field (1965) and Defouw (1970), has been studied in detail analytically by N10, but only for 
< I | unmagnetized media (see, e.g., D'Ercole & Ciotti 1998 and Li & Bryan 2012, for a numerical approach). 

The rotational properties of the hot gas of galaxies and galaxy clusters are almost unconstrained observationally, because 
of the relatively poor energy resolution of current X-ray instruments. There are theoretical reasons to expect that rotation 
is especially important for the coronae of disc galaxies (Marinacci et al. 2011), but some models predict significant rotation 
also for the hot atmospheres of elliptical galaxies (e.g. Brighenti et al. 2009) and galaxy clusters (e.g. Lau, Kravtsov, & Nagai 
2009). In any case, even if not dynamically dominant, rotation influences the thermal-stability properties of the fluid (N10), 
as it happens for a subthermal magnetic field (Balbus & Reynolds 2010). There is observational evidence of the presence of 
magnetic fields in both the intracluster medium (Carilli & Taylor 2002) and the gaseous halos of galaxies (Jansson & Farrar 
2012), though the field geometry is still poorly constrained. Notwithstanding these uncertainties on the detailed kinematical 
and magnetic properties of these gaseous systems, it is clear that allowing for the presence of both rotation and magnetic 
fields would represent an important step forward in the attempt to understand the thermal stability of the hot atmospheres 
of galaxies and galaxy clusters. 

Here we focus on the problem of the linear stability of a rotating, weakly magnetized plasma against axisymmetric 
perturbations, accounting for the effects of stratification (i.e. the fluid is in equilibrium in an external gravitational field), 
differential rotation, radiative cooling and anisotropic thermal conduction. Our analysis can be considered a generalization of 
previous studies on the linear stability of magnetized fluids. For instance, the non-rotating case has been studied both in the 
absence (Quataert 2008; Kunz 2011) and in the presence (Balbus & Reynolds 2008, 2010; Latter & Kunz 2012) of radiative 
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losses. The linear stability of a rotating magnetized plasma, in the absence of radiative cooling, has been widely studied in the 
astrophysical literature, with specific applications to accretion discs (e.g. Balbus & Hawley 1991, 1992; Urpin & Brandenburg 
1998; Kim & Ostriker 2000; Brandenburg & Dintrans 2006; Islam 2012; Salhi et al. 2012) and rotating stars (e.g. Fricke 1969; 
Pitts & Tayler 1985; Menou, Balbus, & Spruit 2004; Masada 2011). The stability-analysis techniques developed in these works, 
as well as in studies of rotating unmagnetized media (see N10, and references therein), can be applied to the more general 
problem addressed here. 

It is worth spending a few words on the structure of the magnetic field, which, as already mentioned, is poorly constrained 
observationally. Throughout this work we will assume that, at least from the point of view of the perturbation, the magnetic 
field is ordered. This does not necessarily imply that the magnetic field is ordered over scales comparable to the size of the 
system: the field may as well be tangled, but with coherence length substantially larger than the size of the perturbation. 
If, instead, the coherence length of the field is comparable to or smaller than the size of the perturbation, the system is 
better described by an unmagnetized model (such as that of N10) in which thermal conduction is suppressed by some factor, 
accounting for the effect of the tangled magnetic field (see, e.g., Binney et al. 2009, for a discussion). 

The paper is organized as follows. In Section 2 we describe the plasma model and we derive the general dispersion relation 
for linear axisymmetric perturbations. Stability criteria for previously studied limiting cases are obtained in Section 3. In 
Section 4 we present our new stability criteria. In Section 5 we summarize our main results and discuss the implications for 
the hot atmospheres of galaxies and galaxy clusters. 



2 LINEAR-PERTURBATION ANALYSIS 
2.1 Governing equations 

A stratified, rotating, magnetized atmosphere in the presence of thermal conduction and radiative cooling is governed by the 
following magnetohydrodynamics equations: 

g + V-(pv) = 0, (1) 



Of 



£ + (v.V)v 



P+^) -PV$+^(B.V)B, (2) 



2£ - V x (v x B) = 0, (3) 



dt 
P 

7-1 



5 + "-* 



Hpp'I = -v • Q - pC, (4) 



supplemented by the condition that the magnetic field B is solenoidal (V ■ B = 0). Here p, p, T and v are, respectively, the 
density, pressure, temperature and velocity fields of the fluid, B = |B|, $ is the external gravitational potential (we neglect 
self-gravity), 7 is the adiabatic index, Q is the conductive heat flux, and £ = C(T,p) is the radiative energy loss per unit 
mass of fluid. In a dilute magnetized plasma heat is significantly transported by electrons only along magnetic force lines (see 
Braginskii 1965). In other words, in the presence of a magnetic field the thermal conduction is anisotropic, so the conductive 
heat flux is given by 

XB(B-V)T 

Q = ^ > ( 5 ) 

where \ is the Spitzer electron conductivity that can be expressed as 

X = *T 5/2 , (6) 

with re ~ 1.84 x 10~ 5 (lnA)~ ergs 1 cm 1 K 7//2 , and In A is the Coulomb logarithm (Spitzer 1962). In the following we 
neglect the weak temperature and density dependence of InA, assuming that k is a constant, so \ — x(T) oc y 5 / 2 . In the 
unmagnetized case thermal conduction is isotropic, i.e. the heat flux is 

Q = -reT 5/2 VT. (7) 

Thermal conduction in the presence of a magnetic field is therefore reduced in directions that are not parallel to the field 
lines, and it is null in the direction orthogonal to the field. It follows that the effect of a tangled magnetic field is, in general, 
a suppression of heat conduction: in a first approximation the presence of a tangled magnetic field can be modeled as an 
unmagnetized medium with a scalar thermal conduction (equation 7), but with re reduced by some factor (Binney & Cowie 
1981). However, if the coherence length of the tangled field is substantially larger than the size of the perturbation, the 
anisotropy of the conductivity becomes crucial for the stability properties of the magnetized medium (Balbus 2001; Quataert 
2008). In the following we will focus on the latter case and we will consider the anisotropic heat flux as given in equation (5). 
In cylindrical coordinates (R, <j>, z), neglecting all derivatives with respect to cj> (because we will consider only axisymmetric 
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unperturbed fields and disturbances), the governing equations (1-4) read 

dp 1 ORpvR dpvz _ 
~di R dR ~d~T ~ ' 



dv R Ovr 8vr 
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dB * _ R dV 4> j_ R 9v <t> , 



vrv^, 



idp_ 

pdR 
1 

i? _ _ 4^rp 
1 <9p 9$ 

p dz dz 
dB R 

- Vr- 



a$ 1 



Bl 
R 



dB^ dB$ B R B^ 
dR oz R 



dB 
dt 

P ( d 

1 - 



dv z 
dz 



t +VR M + v ^ 



d_ 

dz 



OR 

B^VR 

R 

8B Z 
OR ' 

ln(pp ' 



■ v R 



vr 



1 

47rp 
dB R 
dz 
dBj, _ 
OR 
dB z 
dz 

RdR 



dB z dB z 

Br ^r +b >^7 



(V-v)Bk, 

dB$ _ v^Br 
R 



1 dB 2 

8np dz 



dz 

(V-v)B„ 

X(T)B R (B ■ V)T 



V • v)B d 



1 

8^p~dR 



R- 



+ 



dz 



X (T)B Z (B- V)T 



B 2 
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(8) 
(0) 
(10) 

(11) 
(12) 
(13) 
(14) 
(15) 



where we have used V • B = in writing the three components (equations 12-14) of the induction equation (3). 



2.2 The unperturbed plasma 

The unperturbed system is described by time- independent axisymmetric pressure po, density po, temperature To, velocity 
vo = (vor, uo<;>, vo z ) and magnetic field Bo — (Bor, Bq^, Bq z ) satisfying equations (8-15) with vanishing partial derivatives with 
respect to t, under the assumption that the plasma is weakly magnetized, in the sense that the parameter /3 = 8irpo/Bo 3> 1, 
where Bo = |Bo|: in other words, the magnetic field is subthermal and dynamically unimportant. Formally, such a steady-state 
configuration requires that, in the unperturbed system, cooling is perfectly balanced by heat conduction, which appears like 
an artificial and unrealistic assumption. However, the steady-state solution can be interpreted more broadly as describing a 
"quasi-stationary" state, which does not evolve significantly over the timescales of interest and is close to hydrostatic and 
thermal equilibrium. So, globally, the cooling time of the system is assumed to be much longer than the dynamical time, 
which is clearly the case for the hot atmospheres of galaxies and galaxy clusters. 

The unperturbed fluid is allowed to rotate differentially with angular velocity tt(R, z) = vo<f,(R, z)/R depending on both R 
and z. Without loss of generality we choose our azimuthal coordinate <f> so that Q > 0. We assume that there is no meridional 
circulation in the background fluid. In principle the velocity components vqr and vq z could be non-null, even in the absence 
of meridional circulation, because a time-independent subsonic inflow of gas can occur if cooling is not perfectly balanced 
by thermal conduction (see, e.g., N10). However, for simplicity, in the present investigation we limit ourselves to the case 
vor = Voz = 0. It must be noted that the Poincare-Wavre theorem (Tassoul 1978), which holds for unmagnetized fluids, 
applies also to our model of magnetized plasma, as long as the magnetic field is dynamically unimportant in the unperturbed 
configuration (/3 S> 1). Therefore, the fluid is baroclinic [i.e. p = p(p, T)] in the general case in which $1 = Q(R, z), and is 
barotropic [i.e. p — p(p)] only in the particular case in which the angular velocity is constant on cylinders [fi = il(R)]. 

A useful relation among unperturbed quantities is the vorticity equation, which is derived from the momentum equations 
and, under the above hypotheses, can be written as 



R 



dn 2 



fdTodpo 



dT dpo 
dR dz 



(16) 



dz pqTq \ dz dR 

As we have assumed that there are no motions in the meridional plane (vor = vq z — 0) , in the hypothesis that all components 
of the background magnetic fields are time-independent the induction equation implies that the unperturbed system satisfies 
Ferraro's isorotation law (Ferraro 1937) 

B • Vfi = 0, (17) 

i.e. the angular velocity is constant along field lines (see equation 13). It has been noted (Balbus & Hawley 1991) that 
condition (17) might be too restrictive when the magnetic field is weak. If Bo • Vf2 7^ the azimuthal component Bq^ of 
the magnetic field varies secularly (see equation 13), but, provided that the magnetic field remains subthermal, the time 
dependence of Bo,/, does not imply the time dependence of any of the other unperturbed fields. So it is possible to consider a 
more general unperturbed weak-field configuration, in which condition (17) is not satisfied, B04, depends on time, while all the 
other unperturbed quantities (including Bor and Bq z ) are time-independent: the results of the perturbation analysis are still 
valid in this case, provided the dispersion relation does not depend on Bo$. It turns out that this is not the case for the most 
general dispersion relation derived in the present work (see Section 2.3), so we will limit ourselves to the case Bo ■ Vfi = 0. 
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2.3 Dispersion relation for axisymmetric perturbations 

We describe here the linear-perturbation analysis of the weakly magnetized, stratified, dissipative fluid governed by the 
equations reported in Section 2.1, assuming a background configuration as described in Section 2.2. In practice, we linearize the 
system (8)-(15) by using axisymmetric Eulerian perturbations of the form fo + f e ^ iu ' t + lk R R + lk ^ z ^ wnere fo is the unperturbed 
quantity, |/| <C |/o|, <^> is the perturbation frequency, and kn and k z are, respectively, the radial and vertical components 
of the perturbation wavevector. The linear analysis is intended to be taken locally in the plasma, in the sense that the 
perturbation wavelength is much shorter than the characteristic scalelengths of the unperturbed system. We further assume 
that the perturbation frequency is much lower than the sound-wave frequency, so we work in the Boussinesq approximation 
to exclude a priori the (stable) modes describing the sound waves. In linearizing the thermal conduction term, it is useful to 
note that the linear perturbation of the heat flux, in the Wentzel-Kramers-Brillouin (WKB) approximation, is 

Q = -x(2b) [b (b ■ V)T + b (b ■ V)T + b(b ■ V)T - 2b (b ■ b)(b ■ V)T ] , (18) 

where we have introduced the dimensionless vectors b = B/Bo and bo = Bo /-Bo- From the above expression, writing explicitly 
the axisymmetric perturbation, we obtain the thermal-conduction term of the linearized energy equation: 



V Q = x(T )To 



^2 T 



- (k ■ b ) V + i(k ■ b ) ( V In To ■ b) + i( V In To ■ bo) (k • b) - 2i(b ■ k) ( V In T • b ) (b • b) 
Jo 



(19) 



where k = (kn,0,k z ), because k^ = for axisymmetric disturbances. The first term in the right-hand side of equation (19) 
is similar to what is obtained in the case of unmagnetized heat flux (equation 7), but accounts for the relative orientation 
of the displacement and the unperturbed magnetic field. The second term leads to the so-called magnetothermal instability 
(MTI; Balbus 2001). The third and fourth terms lead to the so-called heat-flux driven buoyant instability (HBI; Quataert 
2008): we recall that these latter two terms are null if one makes the assumption of isothermal unperturbed magnetic field 
lines (VTo 'bo = 0), which implies that there is no heat flux in the unperturbed medium. In the present framework, in which 
the gas is allowed to cool, we are interested in the more general case VTo ■ bo 7^ 0, the underlying assumption being that the 
timescale of the background heat flux is long as compared to the local dynamical time. It must be noted that, even when 
radiative cooling is negligible, the assumption of isothermal field lines is not necessary: formally, the unperturbed system can 
be in steady state also when VTo • bo 7^ 0, provided that the divergence of the heat flux vanishes (see Quataert 2008, for a 
discussion). 

Linearizing the system of partial differential equations (8-15), assuming axisymmetric perturbations in the WKB approx- 
imation, we get the following system of linear equations: 

ik R v R + ik z v z = 0, (20) 
-iujvRpo + ik R p - 2Q,v 4> p - c%A pR p - lkzB ° z b r + ^-^(Bo^Bd, + B 0z B z ) = 0, (21) 

47T 47T 

-iuvspo + (0 + Qr)vrPo + tt z v z p - '(k B °) g _ q 5 (22) 

47T 

-\ujv z p + \k z p - c%A vz p - lkRBoR b z + ^-(Bo^Bcf, + BqrBr) = 0, (23) 

47T 47T 

-icuB R - i(k ■ B )vr = 0, (24) 

-iuBtf, - i(k ■ B o )«0 - (Qr - n)B R - Q, Z B Z = 0, (25) 

-\ujB z - i(k ■ B )w* = 0, (26) 

- vyu>— +vr(A p r - jA p r) + v z (A pz — jApz) 

7 L po 

= -oJd^r + i l-^) X(7 ° )r ° [(k ■ b ) V In To + (VlnT. ■ b )k - 2(k ■ b )(VlnT • bo)b ] ■ b, (27) 
io V 7 / Po 

where c§ = po/po is the isothermal sound speed squared; the quantities A p r = (dpo/dR)/po [A p r = (dpo/dR) / po] and 
A pz = (dpo/dz)/po [A pz = (dpo/dz) / po] are the inverse of the pressure (density) scale-length and scale-height, respectively. 
In analogy with N10 we have defined 

Joj c +cj th if Bo = 0, 

LU d = < [AS) 

[ojc.a + wth if Bo 7^ 0, 

where 

k 2 X (T )To 



7 / Po 
is the isotropic thermal-conduction frequency, 

_ (k-bo) 



A- 2 



(29) 



(30) 
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is the anisotropic thermal-conduction frequency, and 

wth = - (— - ) — [C(po,T ) + p £ p (p ,T ) -T Ct(po,T )] (31) 
V 7 / Po 

is the thermal-instability frequency with C p = dC/dp and Lt = dC/dT; Q.r = d(£lR)/dR and fi z = d(Q,R)/dz are two 
frequencies associated with the angular velocity gradient. In terms of the defined quantities, the assumption of short- wavelength 
perturbations gives \kii\, \k z \ 3> |A p .r|, |A pz |, |A pz |, |A p .r|, and f2 2 ,f2fj,f2 z <C c§fc 2 , while the assumption of low-frequency 
perturbations gives to 2 <C c 2 k 2 . As implementations of the Boussinesq approximation we neglected the term —iwp/po in the 
mass-conservation equation and the term — icjp/po in the energy equation, and we assumed poT ~ —Top (see N10). We recall 
that the magnetic field is considered weak (/3 3> 1): in particular we assumed j3 of the order of {k/\A pR \) 2 . 

The system of linear equations (20-27) can be reduced to the following 5-th order dispersion relation for n = — iui: 

rf + uidiT, 4 + [ojbv + ^rot + 2oj a ] n 3 + [(ui 2 ot + 2oj A )cjd + w A cj c ,mag] n 2 + 



( ^A + Wbv + Wr ot - 40 2 ^| + ui 2 !4 , ) n + W A 



kt 2 \ 2 IV 2 . 2 ,^2^\ 2 „ (32) 



^A + ^r 2 ot - 4 ^ 2 7 c t ) ^ d + W A^c,mag 



0. 



This is the most general dispersion relation derived in the present work, describing the evolution of axisymmetric perturbations 
in a gravitationally stratified, rotating, plasma subject to radiative cooling and thermal conduction, in the presence of a weak 
magnetic field of arbitrary geometry. Here, as in N10, we have introduced the quantity 

w «t = -§^V(R 4 n 2 ), (33) 

which is the square of the frequency associated with differential rotation, and the Brunt- Vaisala (buoyancy) frequency cjbv, 
defined by 

2 _ klVp 

Wbv = -77 ©so, (34) 

k poj 

where so = InpoPo" 7 ^ s ^ ne unperturbed specific entropy and 

k R d d 

V =YlYz^dR (35) 
is a differential operator which takes derivatives along surfaces of constant wave phase (Balbus 1995). In addition we have 
introduced the Alfven frequency uia = k ■ va, where va = Bo / y/ iirpo is the Alfven velocity, and the following frequencies 
related to the anisotropy of the thermal conduction due to the presence of a magnetic field: 



47rpo k 2 Vlnpo 



B 2 k 2 k 2 



Pin T - 2 (VlnTo ■ bo) ( ^-b 0z - b 0R 



k z 

which is in general non-null for any geometry of the unperturbed magnetic field, and w c ,^, defined by 



(36) 



2 _ n ^Pok z V\np 
u c ,<t> = WcJ2_ g2~^2 — £2 — (VlnTo • bo)6o0, (37) 



which can be non- vanishing only if the magnetic field has a non- vanishing azimuthal component (bo^ 0). We note that the 
dispersion relation (32) depends on Bq^ through the quantities 6o</>, boa and &o z , so we will restrict our analysis to unperturbed 
configurations with isorotational magnetic field (Bo ■ Vfi = 0), for which Bq^ is time- independent (see discussion at the end 
of Section 2.2). 

We recall that, in terms of the quantity n appearing in the dispersion relation (32), the perturbation evolves with a time 
dependence f(t) cx e nt , where in general n £ C. Therefore, stable modes are those with Re(n) < and unstable modes those 
with Re(n) > 0. Among the unstable modes it is useful to distinguish between purely unstable modes, having Im(n) = 0, in 
which the perturbation grows monotonically, and overstable modes, having Im(n) 7^ 0, in which the perturbation oscillates 
with exponentially growing amplitude. We also remind the reader that, while ui\ > for all k, the quantities u) 2 ot , wlv an d 
Lj 2 j, can be either positive or negative, also depending on k. A list of definitions of some relevant quantities used in the present 
work is given for reference in Table 1. 



3 STABILITY CRITERIA IN LIMITING CASES: COMPARISON WITH PREVIOUS STUDIES 

Before investigating the stability conditions for the dispersion relation (32) in its general form, it is convenient to discuss some 
simpler cases that have been already studied in the astrophysical literature. In particular we analyze here the limiting cases 
of the dispersion relation (32) in which at least one among oj a , cj t h and Q, is zero. In all cases we account for the effect of 
thermal conduction on the stability of the plasma. This preliminary analysis is propaedeutic to Section 4, where we present 
new stability criteria, because it helps simplify the derivation and the interpretation of the results of the more general cases 
discussed there. In this Section, as well as in Section 4, for a given dispersion relation we obtain necessary and sufficient 
conditions for stability following the approach described in Appendix A, which is based on the Routh-Hurwitz theorem for 
the stability of polynomials. 
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Table 1. List of definitions. 



a _ d In po . _ d In po . _ ainpp . _ 8 In p n 

^Pfl — 3fl ' V — 8z ' fl pB — 8fl ' P z — dz~' 

c « = fe) ' ^^fe^-^J- n = -iw, S0 = lnpopo 7 , 

lg = - Wc 4||a||^£o [©i nTo _ 2(VlnT ■ b ) (^6 , - &oh)] , ^ = ^^^I^ 1 (VlnTo • bo) 



"rot 



\ - 'J C 

= -|l^»(J? 4 n 2 ), wa, = - (2=1) ^ [£(po,T ) + p £p(P0,T ) - T £ t (po,T )], 
swth+WcafifBo^O), u; d = u; th + u; c (iffl = 0), n H = ^J., fl, = a = 



3.1 Unmagnetized, rotating, radiatively cooling plasma (SI 7^ 0, o;a = 0, Wth 7^ 0) 

When the Alfven frequency uja = the dispersion relation (32) reduces to 

n 3 + n 2 uj d + (oj| v + u 2 ot )n + a; 2 ot oj d = 0, (38) 

which was obtained by N10 for an unmagnetized rotating plasma. We note that the condition uja = is achieved not only if the 
fluid is unmagnetized, but whenever the projection of the magnetic field onto the wavevector is zero. For example, the dispersion 
relation (38) is obtained also for a system subject to an azimuthal magnetic field (Bo$ 7^ 0, Bqr = Bo~ — 0), perturbed with 
linear axisymmetric disturbances. In this azimuthal magnetic field configuration thermal conduction is ineffective (uid = Wth), 
in contrast with the unmagnetized case in which (w d = ui c + Wth), but in both cases the dispersion relation is formally 
identical to the most general dispersion relation found in N10 for axisymmetric perturbations. Performing a stability analysis 
of equation (38) as described in Appendix A, we obtain the necessary and sufficient stability criterion 

OJd > 0, CJrot > 0, OJBV > 0, (39) 

as found in N10 (see figure 1 in that paper). For the above criterion to be satisfied for all wavevectors the following conditions 
must hold: 

^ th >0, Vp o -V So <0, f^ = 0, W£p.>0. (40) 

oz oR 

So, we have stability only for barotropic fluids [SI = S1(R)[ satisfying the Field (aj tn > 0), Schwarzschild (dso/dpo < 0) and 
Raylcigh d(SlR 2 )/dR > criteria. When these conditions are not satisfied we can have either monotonically growing instability 
or overstability (see N10). 

3.2 Magnetized, non-rotating plasma without radiative cooling (SI = 0, u>a 7^ 0, Wth = 0) 

Here the fluid is magnetized but it does not rotate and there is no cooling. This is the case considered in Cartesian coordinates 
by Quataert (2008). We still work in cylindrical coordinates, but, without loss of generality, we assume Bq^ = (so b^ R + &oz = 
1). The dispersion relation (32) becomes 

(n 2 + lu a ) [n 3 + uic^n 2 + (u>\ + w| v ) n + ui A (w c , a + Wc.mag)] = 0. (41) 

The zeroes of the quadratic factor are n = ±iuA: the two oscillatory solutions describing the Alfven waves. Applying the 
stability analysis (see Appendix A) to the cubic factor we first get the following two necessary conditions for stability: 

A 2 > oj c , a a;lv — a; Cima gCJA > 0, (42) 

A3 > wl(w c ,aWBV — Wc.magWA) (Wc,mag + W c ,a) > 0, (43) 

because Ai = aj c , a > for all k. The condition on A2 can be rearranged as 1 



£>Po 



-Vs Q + 27 In To — 2 (V In T ■ b ) f^6o. - b R 
7 \kz 



> 0. (44) 



Using the identity T>sq — 7©lnTo — (7 — l)2?lnpo, it is possible to write the condition (44) as 

1 - (©po) 2 - 2£>p (VlnTo • b ) ( ^b 0z - b R j > 0. (45) 



7 po \k 



1 It is useful to notice that i^ A aj c>mag = — (x! Cia ^§- £o_X> lnpo ^DlnTo — 2 (VlnTo ■ bo) ^^-boz — &o 
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We note that in the absence of rotation the fluid is barotropic, so surfaces of constant density, pressure and temperature 
coincide (see Section 2.2), and, without loss of generality, we can assume po — po(z) and To = Tb(z), so that T> = xd/dz, 
where x = kn/k z . Requiring for stability the validity of the inequality (45) Vx, we get the following two conditions for stability: 

dT dp n dpo/7-ldlnpo ,2 dlnT \ 

bonbon— — =0, — - 2b 0z — > 0. (46) 

dz dz dz \ 7 dz dz J 

When A 2 > the condition A 3 > becomes cj c ,a H~ ^cmag ^ 0, which leads to the following additional condition: 

-(2b 2 2 -l) >0. (47) 



dTb dp 2 



dz dz 

Let us now consider the singular cases (Ai = or A2 = or A3 = 0): it can be shown that the only additional condition 
implied by these cases is the Schwarzschild stability criterion oj| v > 0, which, in the present case, gives 

^£^< i.e. dpo/ 7-ldlnpo dlnTtA 
dz dz dz \ 7 dz dz J 

Summarizing, the necessary and sufficient criterion for stability is 

b 0zb0R ¥°¥° = 0, %° (l=±****L _ 2&^) > 0, (49) 
dz dz dz \ 7 dz dz J 

^£dP°(2&g, -!)><), d^/ 7 -ldlnp _dlnTo\ Q 
dz dz dz \ 7 dz dz J 

When the unperturbed field lines are isothermal (i>oz = 0, boR = 1), the stability criterion is simply 

^ < (for stability, if b 0z = 0), (51) 
dpo 

where we have used To = To(po)> because the fluid is barotropic (in other words the temperature gradient must be opposite 
to the pressure gradient). If this condition is not satisfied we have the MTI (Balbus 2001). When the unperturbed field lines 
are not isothermal (bo z 7^ 0), the stability criterion is 

b R = 0, < < (for stability, if b 0z ± 0), (52) 

a In po 2 7 

so we have stability only when temperature gradient, pressure gradient and magnetic field are parallel, and temperature 
increases for increasing pressure, but with logarithmic slope smaller than (7 — l)/2 r y. If &oz 7^ 0, but the above stability 
conditions are not satisfied, we have the HBI (Quataert 2008). 



3.3 Magnetized, non-rotating, radiatively cooling plasma (S7 = 0, uia 7^ 0, Wth 7^ 0) 



In this case the plasma is magnetized and not-rotating, but, differently from Section 3.2, we consider the cooling term in the 
energy equation (see Balbus & Reynolds 2008, 2010). As in Section 3.2, the fluid is barotropic and, without loss of generality, 
we can consider Bq^ — (i.e. b^ R + b Q Z = 1). The dispersion relation is 



(n 2 + Wa) [n 3 + uidn 2 + (ui\ + wlv) n + w\ (ojd + ^ c ,ma g )] = 0. 
Applying the analysis of Appendix A to the cubic factor, we first find the necessary conditions for stability 

Ai > ^> u d > 0, 



A 2 > 
A :i > — ■ - 2 ' 2 



2 



■ ^c.mag^A ^ 0, 



(53) 

(54) 
(55) 
(56) 



> (^(^d^BV — W Cil]la gUA)(Wc,mag + Wd) > 0. 

Imposing the validity of (54) for every wavevector k, we find u t h > 0, i.e. Field (1965) stability criterion. Condition (55) can 
be rewritten as 



/ . \ 2 I k z p 

(W c ,a +LJth)^BV + W c , sT t L>mpo 

' k po 



Z?lnr -2(VlnT -bo) 



-h} Z 



> o, 



which is valid for every wavevector k when 



2 k z p 
k 2 po 



^th^BV > 0, 



PlnT -2(VlnTo -bo) 



k R 



boR 



> 0. 



(57) 

(58) 
(59) 



As in Section 3.2, we can assume po = po(z) and To = Tq(z), so T> — xd/dz, where x = kR,/k z - Given that u t h > 0, 
condition (58) reduces to the classical Schwarzschild criterion u) BV > 0, which in this case can be written as condition (48). 
For the inequality (59) to be satisfied for all k we must have 



dp dT 

— ; j — bo z boR = 

dz dz 



and 



dpo 
dz 



7 " 



ldlnpo din Tb ,2 

* : O 0; 



d; 



dz 



> 0. 



(60) 
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When A2 > condition (56) can be simplified in + w Cimag > 0, which gives 



dp dT . 2 v 2 . o dp dT , . 

-3 3— (l-2b 0z )a: + 2— — b 0z b 0R x < 0. (61) 

ds dz dz dz 

For this to be true for all x we have the additional condition 

*S*p-aSD<a («> 

It can be shown that the singular cases (Ai = 0, A2 = and A3 = 0) do not lead to additional stability criteria. Summarizing, 
in the present case the necessary and sufficient stability criterion is 

dp f -y-1 dlnpo d In T \ dp dT 

uth > 0, -3— > 0, — — feozOofl = 0, (63) 

dz \ 7 dz dz ) dz dz 

dpo^-ldlnp 2 dta3b^> Q) ^£d7o (1 _ 26L) < . (64) 



dz \ 7 dz dz J dz dz 

When the unperturbed field lines are isothermal (t>oz = 0) we get 



wth>0, t^<0 (for stability, if b 0z =0); (65) 
dpo 

when the unperturbed field lines are not isothermal (bo* 7^ 0), we get: 

wth>0, 6 fl = 0, 0<^-^<^— - (for stability, if b 0z / 0), (66) 

dlnpo 27 

in agreement with Balbus & Reynolds (2010). So, formally, we can have stability for specific magnetic field orientations. 
However, as well known, the Field criterion (oj t h > 0) is typically not satisfied in astrophysical plasma. When the Field 
criterion is violated, while the other stability criteria are met, either overstability or monotonically growing instability occurs 
(see Balbus & Reynolds 2010). 



4 NEW STABILITY CRITERIA FOR A WEAKLY MAGNETIZED, ROTATING, STRATIFIED 
PLASMA 

We are now going to analyze the dispersion relation (32), describing the evolution of linear perturbations in a stratified, weakly 
magnetized, rotating plasma subject to radiative cooling and thermal conduction. Before addressing the most general case, 
which we will discuss in Section 4.4, it is convenient to consider separately cases in which some terms appearing in equation (32) 
vanish. In all these cases, which, as far as we are aware, have not yet been treated in the astrophysical literature, we account 
for the fact that the medium is magnetized (uja 7^ 0), rotating (fi 7^ 0) and subject to thermal conduction along the field 
lines, and we always allow for the unperturbed magnetic field lines to be non-isothermal (VTb - bo 7^ 0), but we restrict our 
analysis to systems with isorotational unperturbed magnetic field lines (bo • V$l = 0). 



4.1 Plasma in the absence of radiative cooling (wth = 0) with meridional magnetic field (£>o<£ = 0) 



Here we derive stability criteria for a rotating, magnetized plasma in the absence of radiative cooling. For simplicity we assume 
here that the magnetic field has no azimuthal component, i.e. _Bo</> = (therefore, b^ R + bf Sz = 1). Neglecting radiative losses 
we get the dispersion relation 

n 5 + cj c , a n 4 + (wlv + ^rot + 2o;i) n 3 + [(wrot + 2a;i)a; c>a + ajiaj c , mag ] n 2 + 



. 2 / 2 . 2 , 2 . n 2 n-z \ . 2 



0. 



Analyzing this dispersion relation as in Appendix A we first derive the following necessary conditions for stability: 



A 2 > 



A 5 > 



> 0, 



^c,a ( + <*>rot — 4f2 -| J + l^A^c.mag > 0, 



(67) 



(68) 



(69) 



because Ai and A4 are always nonnegative (A4 = only if A2 = 0), and A3 > if A2 > and A5 > 0, being A3 oc 
A 2 [wc,a (oj% + Wf 0t ) + w A w c ,mag] • The condition (68) gives 



7-1 



7 



V\np - 2(VlnT • b ) 



bv>z — boR 



> 0. 



For this to be satisfied for all wavevectors we should have 
[7-I 



Vp 



7 



-Vlnpo-2(VlnTo-b )b 



> 



and (VT • b ) (Vp x b ) < 0. 



(70) 
(71) 



Stability of a rotating plasma 9 



Of course the latter condition can be satisfied only with the equality, so we can rewrite it as 

(VT -bo)(Vp x b ) = 0. 
Let us now consider the condition (69): the criterion to have A5 > for all wavevectors is 



-Vpo • VlnTo +p 



dlnR 



+ 2(VlnTo-b )(Vpo-bo)>0, 



dpo 

dz 



dndtf dT dQ 2 \ f on 2 on 2 



(VT • b ) 2 (Vp x b ) 



> 0, 



(72) 



(73) 



(74) 



RpoTo 

where we have exploited the vorticity equation (16). The singular cases to be analyzed separately are A2 = 0, A3 = 
and A5 = 0. We note that A2 = if either Dpo = [i.e. x — (dpo/dR)/(dpo/dz), if the pressure gradient is non-null] or 
x = —boz/boR- When x = (dpo/dR)/(dpo/dz) the dispersion relation is 



(n + a; Cia ) 



4 , 1 2 n 2 n 2 2/2 2 , n 2 s j 

n + (Wrot + 2oj A )n + uj a [ u TOt +uj a ~ 4S2 — 



= 0. 



(75) 



The mode associated with the linear factor (n = — aj c , a ) is stable because uj c ^ > (conduction damping). Imposing n 2 < in 
the quart ic factor, we get 

dpo ( dQ 2 dpo 9fi 2 dpo " 



dz \ OR dz 



dz dR 



> 0. 



(76) 



When x = —boz/boR the dispersion relation is n 2 + uj 2 ot + ajj| v = 0, so we have stability when D 2 ot + Wbv — 0j which in the 
considered case (a; = —b az /boR) gives 

7 - 1 (b ■ Vp ) 2 



- — (b • Vp ) (b ■ V In T ) + b 0R 

po 



b ■ V(Q. 2 R 4 



R 3 



> 0. 



(77) 
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It can be shown that the singular cases A3 = and A5 = do not introduce additional conditions for stability with respect 
to those obtained above. 

Summarizing, in this case the necessary and sufficient criterion for stability is 

[7-1, 



Vp 



7 



-Vlnpo-2(VlnT -bo)b 



> 0, (VTo ■ b ) (Vpo x b ) = 0, 



-Vp • VlnTo + po 



do. 2 



dpo 
dz 



dlnR 

dTa_Xf_ _ dT dtt 2 
dz dR dR dz 



+ 2(VlnT ■b )(Vp„ -b ) > 0, 



2 (VTo ■ b ) ( boz^z — b R^— 
dR dz 



> 0, 



(78) 

(79) 
(80) 



dpo fdo 2 _dp 1 _xfdpo s . 

dz I dR dz dz dR ' 



7 1 (b ■ Vp ) 2 _ 1_ (bo Vpo) (bo . vinTo) + 46g fl n 2 > 0, (81) 
7 P0P0 Po 



where we have used the isorotation condition bo ■ Vfi = 0. It is important to note that the second of conditions (78) requires 
that, for stability, we must have either VTo ■ bo = (isothermal field lines) or Vpo x bo = (isobaric surfaces orthogonal to 
the field lines), so we now consider these two possibilities. 



4-1.1 Stability criterion when VTo • bo = 

When VTo • bo = 0, in the current hypothesis of isorotation (bo ■ Vfi = 0), the isothermal and isorotational surfaces are 
parallel to the magnetic field lines, so we have fi = fi(Tb) and the dispersion relation simplifies considerably. Analyzing the 
dispersion relation in this limit, we get the following necessary and sufficient stability criterion for a rotating, magnetized 
medium with isothermal unperturbed field lines: 

dp dT 



-Vpo- VlnTo + >0, 
po d m R 



dz dz * °' ( 82 ) 
where we have used the vorticity equation (16). So a necessary condition for stability is that the vertical gradients of temper- 
ature and pressure are opposite (we recall that in this case isothermal and isobaric surfaces do not necessarily coincide). 



4-1-2 Stability criterion when Vpo x bo = 

In this case Vpo x bo = and bo • Vfi = 0, so isobaric surfaces are orthogonal to the field lines and to isorotational surfaces, 
which implies, for instance, that (VTo • bo)(Vpo • bo) = VTo • Vpo- The necessary and sufficient criterion for stability is 

dQ 



7- 



7 



1 |Vlnp | 2 - 2VlnT ■ Vlnp > 0, 



<9fi 2 

Vpo-VlnTo + poTrr-5 >0, 
m it 



dR 



>0, 



(83) 
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Figure 1. Stability and instability domains for rotating, barotropic [fi = Q(R)] fluids with isorotational unperturbed magnetic field 
lines (Bo ■ Vfi = 0, so Bqr = 0), in the space of two parameters: angular-velocity gradient (x axis) and gradient of temperature with 
respect to pressure (y axis). Panel (a) shows the case of a plasma with isothermal unperturbed magnetic field lines (VTo ■ Bo = 0, 
so dTo/dz = 0) in the absence of radiative cooling (tJ t h = 0; see Sections 4.f and 4.2). Panel (b) shows the case of a radiativcly 
cooling, Field-stable (u>th > 0) plasma with isothermal unperturbed magnetic field lines (see Sections 4.3 and 4.4). Panel (c) shows the 
case of a plasma with purely meridional unperturbed magnetic field (Bq^ = 0), with field lines orthogonal to the isothermal surfaces 
(VTq X Bq = 0, so dTg/dR = 0), either in the absence of cooling (aj t h = 0; see Section 4.1) or in the presence of cooling, but with lo^ > 



(see Section 4.3). In panels (a) and (b) dlnf2/dlnii is normalized to (7 ■ 
represent d In To /d In po = a d In Q/d In R. 



1)7" 



where a = 2f2 2 poPo/|Vpo| 2 , so the diagonal lines 



> 0, 1 1^5°L - Vpo ■ V In To + 4p b 2 0R Q 2 > 0, 

7 Po 

(84) 

so a necessary condition for stability is that the angular velocity increases for increasing R. 



dpo 
dz 



dz 8R 



dToXr 
dR dz 



on 2 



dR 



dz 



4-1-3 Stability criteria for a barotropic fluid 

While in the general case of a baroclinic fluid the stability conditions (78-81) are quite involved, in the special case of a 
barotropic fluid the stability criterion becomes considerably simpler and can be easily shown graphically. In the barotropic 
case, the necessary condition (VTo ■ bo) (Vpo x bo) = implies that for a configuration to be stable the magnetic field lines 
must be either parallel or orthogonal to the isothermal surfaces [which are also isobaric because To = To(po)]- Note that, when 
fl — Q(R) the isorotation condition (17) implies that the magnetic field is vertical, so here we must assume bon = 0. 

When the unperturbed field lines are isothermal (VTo • bo = 0, which in the current hypotheses implies dTo/dz = 
dpo/dz — d&/dz = 0) the stability criterion reduces to 

din To dlnfi 2n 2 p oPo 

— < a— - — — , where a = —= — r— (85) 

dlnpo ~ dlniT |Vp | 2 

is a positive dimensionless factor. The corresponding domains of stability and instability are visualized in Fig. 1(a). Note 
that, in this special case, there can be stable configurations with dQ/dR < 0, provided the temperature-pressure gradient is 
strong and negative. In the limit of uniform rotation (dSl/di? = 0) we recover the condition for stability against the MTI (see 
Section 3.2). 

When the unperturbed field lines are parallel to the temperature gradient (VTo x bo = 0, which in the current hypotheses 
implies dT /dR = dpo/dR = and d$/dR = fl 2 R) the stability conditions (78-81) reduce to 

which are visualized in Fig. 1(c). So, a barotropic fluid with magnetic field lines orthogonal to the isothermal surfaces is stable 
if and only if the angular velocity increases outwards and the temperature decreases in the direction of decreasing pressure 
with sufficiently shallow gradient. In other words, the condition of outward increasing angular velocity must be added to the 
condition for stability against the HBI found for non-rotating media (see Section 3.2). We recall that in typical astrophysical 
systems the angular velocity decreases outwards: the fact that the condition dQ/dR. > is violated is at the heart of the 
magnetorotational instability (MRI; Balbus & Hawley 1991, 1992), which is believed to be one of the key mechanisms at work 
in accretion discs. 
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4.2 Plasma in the absence of radiative cooling (wth = 0) with non-vanishing azimuthal magnetic field 
component (B^f, 7^ 0) 

We now discuss the case in which the plasma does not cool, but the unperturbed magnetic field has a non-vanishing azimuthal 
component (Bo</> 7^ 0). A similar analysis was carried out by Balbus (2001), who assumed isothermal unperturbed magnetic 
field lines (VTo ■ bo = 0): we generalize Balbus' result, allowing for non-isothermality of the background field lines. When 
VTo ■ bo 7^ and B04, 7^ 0, but Uth = 0, the dispersion relation is 

n + LOc.a.n* + [w| v + uj 2 ot + 2uj a ] "n 3 + [(^rot + 2u A )o; Ci a + W A w c ,mag] n 2 + 



2 f 2 
+UJ A UJ A 



2 

LOBV 



-4fr 



^2 +^c,4, )n + uj A 



2 2 . n 2*z 



.ma.fi, 



= 0, 



(87) 



which is obtained from equation (32), substituting lu^ with oj c>a (because oj t h = 0), and differs from equation (67) for the 
presence of a term proportional to to c ^ oc (VT • bo) 2 B(^, 7^ 0. For stability, we first impose, following Appendix A, that all 
the Hurwitz determinants are nonnegative. We get: 



A 2 > 



A 3 > 



2 



;^A > 0, 



A 4 > 



2 



Uc,ai^>A + <^c,mag<^>A + W c ,aWrot) — w c> w c,aWA > 0, 
■W?, a wIw C) ^ + (wc.magWA + Wc.aWrot) A 2 tJc,0 + 4fi 



i- 2 
2 K z a 2 



>o, 



(88) 
(89) 

(90) 



A 5 > 



id A + ^r 2 ot - 4f2 2 



A-' 



> o, 



(91) 



(92) 



because Ai > always. Let us consider the condition A4 > 0. This can be written in the following form 

-k 2 f(x) + g(x) > 0, 

where x = kn/k z and 2 f(x) — (oj c>a oj A o; 2 ^/A: 4 ) 2 > for all x. Necessary condition for the above inequality to be satisfied 
for all k (i.e. all x and k) is f(x) = for all x. In the current hypotheses, this implies i>o.R = b 0z = 0, so k • b = for all 
k. However, if b^R = bo* = then uj c ^ = 0, in contrast with our hypothesis, so we must conclude that there is no stable 
configuration with VTo ■ bo 7^ and Bq^ 7^ 0. 

In the case in which the unperturbed field lines are isothermal and Bo^ 7^ 0, u) c>< f, = and the dispersion relation is again 
given by equation (67) with VTb • bo = in a> c ,mag, which is the same as that derived, under the same hypotheses, by Balbus 
(2001). Analyzing this dispersion relation in the same way as in Section 4.1, we get the following stability conditions: 

VTo ■ b = 0, (93) 

~)2 O/TI Q n 2 



-Vp - VlnTo + po 



dlnR 



> 0, 



dp fdTo_Xf_ _ dTp dQ. 2 
dz \ dz dR dR dz 



> 0, 



dpo fdn 2 dp 

dz V dR dz 



dQ, 2 dp 
dz dR 



> 0, 



7 - 1 (b ■ Vp ) 



b ■ V(0. 2 R 4 



R 3 



> 0. 



(94) 



(95) 
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The inequalities (94) coincide with those found by Balbus (2001), while the two conditions (95), which follow from the 
analysis of singular cases of the dispersion relation (67), were not reported before and, as far as we can see, are generally 
independent 3 of the two conditions (94). Specializing to the case of isorotational unperturbed field lines, i.e. imposing Ferraro's 
law (equation 17), the conditions (93-95) reduce to the following necessary and sufficient stability criterion: 



VTo • b = 0, 



-Vp 'VlnTo+po 



^i>0, ^£^£<o, 
d\nR dz dz 



(96) 



where we have used the vorticity equation (16). In the case of a barotropic distribution with isorotational and isothermal field 
(so boR = dpo/dz — 0, but &o</> 7^ 0), the stability criterion is again condition (85), which is represented in Fig. 1(a). 



4.3 Radiatively cooling plasma (cj t h 7^ 0) with meridional magnetic field (B, 



04, 



0) 



We move now to the case in which the plasma cools radiatively (wth 7^ 0): for simplicity we start here analyzing the dispersion 
relation 

n + uj d n 4 + [wlv + oj 2 ot + 2w A ] n 3 + [(uj 2 ot + 2uj 2 a )u) A + W A w c ,mag] n 2 + 



+cj a ( u A + oj bv + w rot - Ail — ) n + oj a 



w A + w, ot - 40 2 ^| ) 



= 0, 



(97) 



2 It can be shown that / depends on the wavevector k only through x = kR/k z , so it is independent of the wavevector modulus k. 

3 The conditions (95) must be considered in addition to the conditions (94) only for baroclinic fluids [Q = tt(R, z)], because when the 
fluid is barotropic [fi = n(i?)] the two inequalities (95) reduce to dQ/dR > 0, which is implied by combining the two inequalities (94). 
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which is obtained from (32) in the hypothesis that the unperturbed magnetic field is meridional (Bo^ = 0), so uj c ^ — 0. 
Performing a stability analysis as described in Appendix A, we get the following necessary conditions for stability: 

Ai > ^> u d > 0, (98) 



A 2 > 



A 5 > 



2 2 
CO A + ^rot 



AD, -| J cj d + oj c , mag a; A > 0, 



(99) 
(100) 



because A4 is always nonnegative and A3 = A2 \uJd (wl + w rot) + ^A^c.mag] > if Ai > 0, A2 > and A5 > 0. Imposing 
the above conditions for all wavevectors leads to the following criteria. The condition on Ai gives wth > 0. The condition on 
A 2 gives 



i.e. ^ = 



(VT • b ) (Vp x b ) = 0, Vp 



7-1 



Vp • Vs < 0, 



Vlnp -2(VlnTo-bo)b 



> 0. 



Exploiting the fact that we must have oj t h > (Field criterion) for the condition on Ai, the condition on A 5 gives 

i9f2 BQ 2 

— >0, -Vpo- VlnTo + po^-^ + 2(VlnT -bo)(Vpo-bo) > 0, 



dpo 
' dz 



dTp on 2 
dz OR 



dT dQ 2 
OR dz 



dQ. 2 dQ 2 
b0 *l)R-- b0R lt 



(VT • b ) 2 (Vp x b ) 2 



> 0. 



(101) 
(102) 

(103) 
(104) 



RpoTo 

The singular cases to be treated separately are Ai = 0, A2 = 0, A3 = 0, A4 = and A5 = 0, but it can be shown that they 
do not lead to additional conditions. In summary, taking into account that one of the conditions (101) states that the fluid 
must be barotropic, the necessary and sufficient stability criterion can be written as 

dQ. 



w th >o, n = n(R), (vr -b )(Vp xb ) 



0, 



Vpo 



-Vlnpo -2(VlnT ■ b ) bo 



> 0, -Vpo- VmTo+po 



■-5— -5 2(VT ■ bo) 60: 

az oz 



> 0, Vp ■ Vs < 0, (105) 
2(VlnTo -b )(Vpo -b ) > 0, (106) 
> 0. (107) 



dR 

dQ 2 
din R 



The above conditions must be supplemented by the isorotation law bo ■ Vfi = 0, which for the relevant barotropic case implies 
boR — (so bo? = 1). For a barotropic fluid To — To(po), so the necessary condition (VTb ■ bo) (Vpo x bo) = means that 
for a configuration to be stable the magnetic field lines must be either parallel or orthogonal to the isothermal (and isobaric) 
surfaces. Let us consider first the case in which the magnetic field lines are isothermal (VTo ■ bo = 0, so dTo/dz = dpo/dz = 0, 
because boR = 0). Performing the same analysis as above, in this limit, we end up with the following necessary and sufficient 
stability criterion: 



cj th >0, Sl = il{R), ^>0, 

ClIX 



dlnTo < 7-I 



din To dlnfi 
< a- 



(when VT ■ b = 0), 



(108) 



dlnpo 7 dlnpo din J? 

where a > is the dimensionless factor introduced in equation (85). These conditions are represented graphically in Fig. 1(b). 
When, instead, the magnetic field lines are orthogonal to the isothermal surfaces (VTo x bo = 0, so OTo/dR — dpo/dR = 0, 
because boR = 0), the stability criterion is 



w t h > 0, 
which is visualized in Fig. 1(c) 



n = n(fl), ^>o, 



0< dTnTo < 7-l 



dlnpo 



(when VT x b = 0), 



(109) 



4.4 Radiatively cooling plasma (cj t h 7^ 0) with non-vanishing azimuthal magnetic field component (Bo^ 7^ 0) 

We now discuss the case in which the unperturbed magnetic field has a non-vanishing azimuthal component (Bq^ 7^ 0) and 
the plasma is allowed to cool. This is the most general case studied in the present work, so the dispersion relation is given by 
equation (32), under the assumption that all terms are non-null (in particular, u} t h 7^ in wa, VTo ■ bo 7^ and So</> 7^ 0, so 
<j-> c .4> 7^ 0). For stability let us first impose, following Appendix A, that all the Hurwitz determinants are nonnegative. We get: 

Ai > ^> u d > 0, (110) 

A 2 > Wd^BV — Wc.magWA > 0, (111) 
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A 3 > 



A 2 (w d ^A + ^c.mag^A + ^d^rot) ~ u\^u\bl\ > 0, 



A 4 > 



2 



-Wd^ia;^ + (wc,magWA + ^d^rot) A 2 Wc,0 + 4f^^§A 2 



> o, 



A 5 > 



cj a A 4 



2 | 2 ,,02 I 

^A + W rot - 4S2 — ] CJ d + Wc,magW A 



> 0. 



(112) 
(113) 

(114) 



As done in Section 4.2, we note that imposing the condition A4 > for all wavevectors, in the considered case in which 
botf, =fc 0, we get boa = &oz = 0. But if boa = &o z = then VTo • bo = and ui c ^ = 0, in contrast with our hypothesis, so we 
must conclude that there is no stable configuration with VTo ■ bo 7^ and Bo$ 7^ 0. 

Let us thus consider the case VTo ■ bo = and Bo^ 7^ 0: now o; c ,0 = 0, so the dispersion relation is equation (97) and the 
corresponding stability conditions are 

<90 „ dQ, 

t)z 



VTo ■ b 







0. 



Wth > 0, 



0, 



> 0, 



Vpo • Vs < 0, 



-Vp ■ VmTo+po- 



dQ, 2 



>0, 



dpo 
' dz 



dR 

dTp dtt 2 
dz dR 



(115) 



dn do, 2 

dR dz 



> 0. 



(116) 



din R 

As in Section 4.3, a necessary condition for stability is f2 = Q(R): note that this is a necessary requirement also in the absence 
of magnetic fields (see Section 3.1), but not in the absence of radiative cooling (see Sections 4.1 and 4.2). Let us focus on the 
case of isorotational unperturbed magnetic field lines: combining the conditions (115-116) with Ferraro's law (equation 17) 
we get that the necessary and sufficient stability criterion is 

n o arm d " > n dlnTp 7 -l dmT dlnO 
0, LJth > 0, it — ll(R), -j-p; > 0, — — — < , - ■ — < a - — — , (117) 



VTo ■ b 



dR dlnpo 7 dlnpo din R 

where a > is the dimensionless factor introduced in equation (85). The corresponding domains of stability and instability 
are shown in Fig. 1(b). 



5 SUMMARY AND CONCLUSIONS 

5.1 Magnetothermal stability criteria for rotating stratified plasmas 

In the attempt to make a step forward in understanding the thermal stability of the hot atmospheres of galaxies and galaxy 
clusters, we have presented new stability criteria for a gravitationally stratified, rotating, radiatively cooling, weakly mag- 
netized plasma in the presence of thermal conduction. We found that for such a medium to be stable against all linear 
axisymmetric disturbances it must obey several restrictive conditions, which are reported in Section 4.3 (for meridional mag- 
netic field; equations 108-109) and Section 4.4 (for magnetic field with non-null azimuthal component; equations 115-116), 
and are represented in Fig. 1. Specifically, in all cases for stability the fluid is required to be Field-stable (uj t h > 0), barotropic, 
with outward increasing angular velocity (dfl/dR > 0). Additional requirements concern the relative orientation of the mag- 
netic field lines and the thermal gradient (which must be either parallel or orthogonal) and the gradient of temperature 
with respect to pressure (which has different effects, depending on the relative orientation of the field lines and the thermal 
gradient). The above conditions can be seen as a combination and generalization of criteria for stability against the MRI 
(Balbus & Hawley 1991), the MTI (Balbus 2001), the HBI (Quataert 2008) and the radiatively driven overstability of Balbus 

6 Reynolds (2010). As well known, even just the two conditions oj t h > and dQ/dR > are typically not satisfied in standard 
astrophysical conditions, so we expect always to find at least one axisymmetric mode that is either monotonically unstable or 
overstable. Therefore, at a formal level, we must conclude that a radiatively cooling, rotating, weakly magnetized astrophysical 
plasma is not expected to be stable against axisymmetric perturbations. Similar conclusions were reached for magnetized (but 
non rotating) fluids by Balbus & Reynolds (2010) and for rotating (but unmagnetized) fluids by N10: in both cases either 
overstabilities or monotonic growing instabilities were expected for typical configurations. In other words, the calculations of 
the present paper have shown that the combination of differential rotation and a weak ordered magnetic field does not have 
a stabilizing effect against thermal axisymmetric perturbations. 

In addition to the above results on thermal stability, we have also generalized previous studies of the stability of magnetized 
media in the absence of radiative cooling: in particular we have extended the study of Balbus (2001), by allowing for the 
presence of non-isothermal unperturbed field lines, and the study of Quataert (2008), by including the effect of differential 
rotation (Sections 4.1 and 4.2). The bottom line of such an analysis is that isothermality of the unperturbed field lines is 
a necessary requirement for stability if the magnetic field has non- vanishing azimuthal component, but, when the magnetic 
field is meridional, stable configurations are formally possible also with (non-isothermal) field lines orthogonal to the isobaric 
surfaces. In any case, the general conclusion is that, even in the absence of cooling, either overstability or monotonic instability 
is expected, because some of the conditions for stability are very unlikely to be satisfied in real systems (e.g., dQ,/dR > in 
the barotropic case). Formally the conditions for stability are more restrictive when cooling is effective: while in the presence 
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of cooling the distribution must necessarily be barotropic in order to be stable, in the absence of cooling it is possible, at 
least in principle, to have baroclinic distributions stable against axisymmetric perturbations. However, these configurations 
are unlikely to be astrophysically relevant, because of the requirement of outward increasing angular velocity. Therefore it 
appears that radiative cooling is not the key factor determining the instability of a rotating magnetized plasma. 

5.2 Implications for the hot atmospheres of galaxies and galaxy clusters 

At a formal level, the above conclusions are the main results of the present work. A quite different question is what these 
results imply from the astrophysical point of view: in particular, what are the implications for the evolution of the hot 
atmospheres of galaxies and galaxy clusters? We recall that the basic underlying question is whether cool clouds can condense 
out of a hot, stratified medium as a consequence of the instability. To answer this question one must necessarily go beyond 
the formal conclusion that linear instability is expected and try to investigate the nature of the instabilities (monotonically 
growing or overstable), their non-linear evolution and the dependence on the properties of the perturbation. Unfortunately, 
we cannot address quantitatively these points based only on the presented calculations, but it is worth discussing them at 
least qualitatively. 

As first pointed out by Malagoli et al. (1987), an important aspect of the problem is the study of the unstable modes, and, 
in particular, determining the nature of the instabilities: while monotonically growing instabilities are expected to naturally 
lead to condensation, overstable disturbances are half the time overdense and half the time underdense, with respect to the 
local background, so they are likely disrupted by turbulence before condensing significantly (Binney et al. 2009; Joung et al. 
2012). The problem of the nature of the instabilities found in the present work could in principle be tackled analytically by 
studying the sign of the real roots of the dispersion relations, but, as mentioned in Appendix A, given the high order of the 
involved polynomials the calculations are prohibitively cumbersome (Liang & Jeffrey 2009, and references therein). Though a 
full exploration of the parameter space is impractical, in selected cases the question of the nature of the thermal instabilities 
can be addressed either by solving numerically the dispersion relations derived in this work or with magnetohydrodynamics 
numerical simulations, which can also be used to gain insight on the non-linear evolution of linearly unstable perturbations. 
For the non-rotating case, an attempt in this direction has been recently made by McCourt et al. (2012): the results of their 
somewhat idealized magnetohydrodynamics simulations suggest that no significant condensation occurs in the relevant regime 
in which the cooling time is much longer than the dynamical time. 

It must be noted that all the stability criteria reported in the present work concern stability against axisymmetric 
perturbations. In principle, configurations that are stable against axisymmetric perturbations could be unstable when more 
general (non-axisymmetric) perturbations are considered. Therefore, from a mathematical point of view, exploring different 
kinds of disturbances could just lead to even more restrictive stability criteria. However, from the astrophysical point of 
view, it is clear that an axisymmetric mode is not necessarily relevant to real disturbances, and a natural question to ask is 
whether the instability occurs also for possibly more relevant non-axisymmetric perturbations. As well known, in the presence 
of differential rotation, it is considerably more complex to study non-axisymmetric than axisymmetric perturbations (e.g. 
Cowling 1951; Balbus & Hawley 1992), and such a task is beyond the purpose of the present paper. However, it is interesting 
to note that, at least in the unmagnetized case, non-axisymmetric perturbations tend to be stabilized by differential rotation, 
while this is not the case for axisymmetric disturbances (N10). 

Throughout the paper we have assumed that the magnetic field is ordered over scales larger than the perturbation size. 
The above-presented calculations show that, if this is the case, the fluid is unlikely to stable, at least against axisymmetric 
disturbances. It is hard to predict the outcome of such instability, but in any case we should expect that also the structure of 
the magnetic field is affected. If turbulent motions are produced by the instability, the (weak) magnetic field is likely tangled 
by turbulence, and the system might end up in a configuration which is better described by an unmagnetized medium with 
suppressed isotropic thermal conductivity (Binney et al. 2009), so that the results of N10 should apply. However, it is also 
possible that, as a consequence of the instability, the fluid rearranges the magnetic field lines in a special ordered configuration 
that tends to contrast the instability, possibly leading to an overstable state (Balbus & Reynolds 2010). In the absence of 
realistic non-linear calculations it appears difficult to decide which of the two above hypotheses is a better description of real 
systems, and of course alternative scenarios are not excluded. Overall, as far as the astrophysical implications are concerned, 
the results of the present work, though not allowing us to draw general conclusions on whether cool gas clouds can condense 
spontaneously out of an equilibrium plasma, strongly suggest that the interplay of rotation and magnetic fields is important 
for the thermal instability, thus encouraging further theoretical and observational investigations of the magnetic and rotational 
properties of the hot atmospheres of galaxies and galaxy clusters. 
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APPENDIX A: ANALYSIS OF THE DISPERSION RELATIONS: STABILITY CRITERIA FOR 
PARAMETRIC POLYNOMIALS 

The dispersion relations obtained from our linear-perturbation analysis are polynomials in the complex variable n = — icj, 
whose coefficients are real and generally depend on the wavevector k. To obtain the linear-stability criterion we need to find 
the conditions under which the dispersion relation is stable for all wavevectors. From a mathematical point of view, this 
problem reduces to studying the stability of families of polynomials of a complex variable with real parametric coefficients. 

We recall here the Routh-Hurwitz theorem, which is a fundamental tool for studying the stability of polynomials of a 
complex variable (e.g. Gantmacher 1959). Let us consider an m-th degree polynomial 

p(z) = aoz m + a\z m ~ l H hOm! ao, Oi, . . . , a m € R, ao=£0, (Al) 

where z £ C and m > 1: p(z) is said to be a Hurwitz polynomial, or simply a stable polynomial, if all its (generally complex) 
roots have negative real part. Let us associate with the polynomial p(z) the following Hurwitz determinants: 



A, EE 



ai a 3 a 5 a 7 
ao a-i a4 a% 
ai a 3 a 5 
ao <Z2 <Z4 



a 2j -i 
a 2 j-2 

a 2 j-4 





for i > m. The Routh-Hurwitz theorem states that when Aj 7^ for j 



(A2) 



1, . . . , m, the polynomial 



where it has been set at 
p(z) is stable if and only if all Hurwitz determinants are strictly positive: 

Ai > 0, A 2 > 0, • • • , A m > 0. (A3) 

A corollary of the theorem is that, under the same hypotheses, a necessary condition for p(z) to be stable is that all its real 
coefficients be strictly positive; i.e. ao, ai, . . . , a m > 0. 

It must be stressed that when at least one of the Hurwitz determinants is null the Routh-Hurwitz theorem does not 
apply. In our application to the dispersion relations the coefficients of the polynomial are parametric and it is not generally 
the case that all the Hurwitz determinants are non-null. Therefore, in order to obtain necessary and sufficient conditions for 
stability, we proceed as follows. We first find the conditions to have 

Ai > 0, A 2 > 0, ■ ■ ■ , A m > 0, (A4) 

for all wavevectors. At this stage these conditions must be considered necessary, but not sufficient, because it is not guaranteed 
that the singular cases (A; = for at least one i) are stable. We then consider separately the dispersion relations obtained 
in the singular cases, thus obtaining additional stability conditions that, combined with those obtained from (A4), give the 
necessary and sufficient stability criterion. 

When the criterion for stability is not met, one is interested in determining whether the unstable modes are monotonically 
growing (z is real and positive) or overstable (z is complex, with positive real part). Mathematically this requires a complete 
root classification of polynomials with real parametric coefficients. This task, which can be easily accomplished for third-order 
polynomials such as those considered in N10, is extremely complex for the higher-order polynomials we deal with in this work 
(e.g. Liang & Jeffrey 2009). Therefore, in the present study we limit ourselves to deriving stability criteria and we do not 
attempt a classification of the nature of the unstable modes. 
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